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We present a new type of transmon split-junction qubit which can be tuned by Meissner screening 
currents in the adjacent superconducting film electrodes. The best detected relaxation time (Ti) was 
of the order of 50 ps and the dephasing time {T 2 ) about 40 ps. The achieved period of oscillation 
with magnetic field was much smaller than in usual SQUID-based transmon qubits, thus a strong 
effective field amplification has been realized. This Meissner qubit allows an efficient coupling to 
superconducting vortices. We present a quantitative analysis of the radiation-free energy relaxation 
in qubits coupled to Abrikosov vortices. The observation of coherent quantum oscillations provides 
strong evidence that vortices can exist in coherent quantum superpositions of different position 
states. According to our suggested model, the wave function collapse is dehned by Caldeira-Leggett 
dissipation associated with viscous motion of the vortex cores. 

PACS numbers: 03.67.Lx, 74.50.-l-r, 85.25.Am 


I. INTRODUCTION 

Several promising architectures for quantum comput¬ 
ers are based on the superconducting qubits based on 
Josephson junctions (JJs). These designs utilize either 
charge^, phase® or fluxed degrees of freedom. These sys¬ 
tems have made tremendous progress in recent years in 
realizing increasingly sophisticated quantum states, mea¬ 
surements and operations with high fidelity®. Supercon¬ 
ducting quits are also attractive technologically because 
they can be naturally integrated into large-scale quantum 
circuits®dO_ However, this main advantage of supercon¬ 
ducting qubits brings a substantial challenge at the same 
time since strong coupling also implies a substantial inter¬ 
action between the qubits and their environment, which 
can break quantum coherence. Understanding the lim¬ 
iting factors of qubit operation is of fundamental and 
practical importance. Previously, a separation of vari¬ 
ous contributing factors to the qubit relaxation and de¬ 
coherence has been achieved^®^®^. Abrikosov vortex is 
one important example of such environment. Dissipation 
caused by vortices has also been studied in a supercon¬ 
ducting resonator^®. Recently, vortices have been shown 
to trap quasiparticles in superconducting resonator and 
qubit, leading to the increase in the quality factor of res¬ 
onator^® and the relaxation time of qubit^^. 

One interesting theoretical possibility of drastically im¬ 
proving quantum coherence in qubits is to couple them to 
Majorana fermions^®. A qubit based on Majorana states 
is expected to exhibit especially long coherence times. 
One approach to create Majorana states is to deposit a 
superconductor onto a topological insulator and to cre¬ 
ate vortices in the superconductor. In this case Majorana 
states can nucleate in the vortex core^®. Thus a study 
of qubits coupled to vortices is needed, in order to deter¬ 
mine whether a qubit coupled to a vortex can preserve 
its quantum coherence and for how long. Here we use 
Meisser qubits to quantify dissipation effects produced 
by single vortices. 


A major advance in the superconducting qubit perfor¬ 
mance became possible after the invention of the trans¬ 
mon qubit^^ . When combined with the three-dimensional 
circuit-quantum electrodynamics (cQED) platform de¬ 
veloped in Ref.^®, the transmon has shown huge improve¬ 
ments of the relaxation time, up to several hundreds mi¬ 
croseconds^ Like the common transmon our device 
involves a capacitance linked by a nonlinear kinetic in¬ 
ductance (JJs). The main difference is that our qubit 
is coupled to the Meissner current and the supercurrents 
generated by vortices. Yet the relaxation time of such de¬ 
vice, designed to probe the environment is rather large, 
namely about 50 ps in the best case. We argue that the 
limiting factor was the Purcell effect, thus the relaxation 
time can be made even longer if necessary. 

Because in our qubit design the Meissner current is 
allowed to flow, partially, into the qubit, a significant 
amplification of the magnetic field effect is demonstrated. 
The qubit transition frequency is periodically modulated 
by the applied magnetic field, but the period is much 
smaller compared to the value estimated by dividing the 
flux quantum by the qubit loop geometric area. 

The Meissner qubit allows a strong coupling to the 
vortices in the leads. We perform a detailed study of the 
radiation-free decoherence effects produced by the vor¬ 
tex cores. It should be stressed that the qubit relaxation 
time may be shortened by the presence of a vortex in the 
superconducting film, due to the Bardeen-Stephan vis¬ 
cous vortex flow. No quantitative study has been done 
so far to test how qubit quantum states relax due to cou¬ 
pling to Abrikosov vortices. Our key finding is that vor¬ 
tices can remain, over many microseconds, in quantum 
superposition states generated because the Lorentz force 
experiences strong quantum fluctuations. The relaxation 
rate added to the qubit by each single vortex was mea¬ 
sured and appears to be surprisingly low, of the order 
of 10 kHz. We propose a semi-quantitative model which 
allows us to estimate this radiation-free relaxation rate 
caused by viscous flow of vortices. Up until now it was 
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FIG. 1. (a) Optical image of the Meissner transmon qubit 
fabricated on a sapphire chip, which is mounted in the cop¬ 
per cavity, (b) A zoomed-in optical image of the qubit. Two 
rectangular pads marked A1 and A2 act as an RF antenna 
and shunt capacitor, (c) Scanning electron microscope (SEM) 
image of the electrodes marked El and E2, and a pair of JJs. 
(d) Schematics of the Meissner qubit. The X, Y and Z de¬ 
note the width, the distance between the electrodes, and the 
distance between two JJs, which are indicated by x symbols. 
The red dot and circular arrow around it in the bottom elec¬ 
trode represent a vortex and vortex current flowing clockwise, 
respectively. 0v is a polar angle defined by two dashed lines 
connecting the vortex and two JJs. The orange rectangular 
loop on the boundary of the bottom electrode indicates the 
Meissner current circulating counterclockwise. 


well established that classical supercurrents can generate 
heat through viscous flow of vortices^®. Now we estab¬ 
lish that quantum superposition currents, such as those 
existing in qubits and characterized by zero expectation 
value, can also generate heat through the same mecha¬ 
nism. Such heat dissipation occurs through the spread 
of the wavefunction of the vortex center followed by a 
collapse of this smeared wave function. 


II. EXPERIMENTAL RESULTS 
A. Qubit frequency modulation 

The design of our devices is shown in Fig. 1, while 
details of the fabrication and measurement techniques 
are described at length in the Appendix A. In brief, the 
qubits have been placed inside a three-dimensional (3D) 
microwave cavity made of Cu (Fig. la). The state of 
the qubit has been determined by measuring the trans¬ 
mission of the 3D cavity. First, we investigated the 
magnetic field dependence of the qubits, anticipating to 
observe periodic SQUID-type oscillations. The trans¬ 
mission versus the magnetic field (“B-field”) varied as 
shown in Fig. 2(a). This plot represents the heterodyne 
voltage, produced by mixing a microwave signal passing 
through the cavity containing the qubit and a reference 
signal. During this measurement the qubit remains in 
its ground state. Yet the cavity input power is chosen 
such that the transmission of the cavity is the most sen¬ 
sitive to the qubit transition from the ground to the ex- 
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FIG. 2. (Golor online) (a) Periodic heterodyne voltage os¬ 
cillation (“HV-oscillation”) as a function of magnetic field. 
Four different colors represent the separate measurement runs 
where the magnetic field swept either round-trip or one-way. 
AB shows the period of the modulation. The arrows marked 
Bi and B 5 indicate the positions of two adjacent sweet spots, 
characterized by dfoi/dB=0. Each sweet spot is equivalent to 
B—0 state. The actual B=0 sweet spot is outside the range of 
the plot. The mark “S” indicates the field at which the peri¬ 
odic signal was shifted to the right, revealing some hysteresis 
caused by vortex entrance to the electrodes. This hysteretic 
behavior can be seen in the second (red) segment of the B- 
field sweep ranging from about 1.45 to 2 Gauss. The mark 
“J” shows the moment when a vortex (or a small group of 
vortices) entered the electrodes during the forward-sweeping 
magnetic field, (b) Comparison of theoretical (AJlff) and ex¬ 
perimental (A“f) effective areas of five samples. The dashed 
line depicts the ideal case of See text for the defi¬ 

nition of these quantities. 


cited state, i.e., the maximum-contrast power was used 
[23] . The probing microwave frequency for this measure¬ 
ment equals the bare cavity frequency. Four segments 
in different color (color online) represent four separate 
measurement runs. The magnetic field was swept round- 
trip (up and down) in the first three segments (black, 
red and green), while it was swept one-way (up) in the 
last segment (blue). The modulation of the transmis¬ 
sion at low iJ-field arises from the change of the onset 
power—the lowest power at which the cavity starts to 
show a sharp increase in transmission, which was termed 
“bright state” (near-unity transmission) in Ref. 23. The 
onset power depends on the difference between the qubit 
transition frequency (foi) and the bare cavity frequency 
(/c). The key point is that /oi is modulated period!- 




























3 


TABLE I. Comparison of the experimental and theoretical 
periods for five samples. AB is the measured period of oscil¬ 
lation. AByz, AByz+cxz, and ABkvz+cxz are the theoretical 
periods calculated using three different effective areas denoted 
by the subscripts. 


Sample X Y 

{pm) (pm) 

Z AB AByz AByz+c 
(^ m) (G) (G) (G) 

;XZ ^-^kYZ + cXZ 

(G) 

N1 

10 

1 

5 

0.38 4 

0.49 

0.30 

N2 

10 

2 

5 

0.3 2 

0.43 

0.29 

N3 

15 

1 

5 

0.43 4 

0.34 

0.22 

N6 

10 

2 

8 

0.16 2.5 

0.30 

0.19 

N7 

10 

2 

8 

0.2 1.25 

0.27 

0.18 


cally by the applied magnetic field, because the qubit 
includes a SQUID-like loop formed by the two JJs and 
the two electrodes, marked in Fig.lc as El and E2. The 
heterodyne voltage Vh is proportional to the microwave 
transmission. Thus, as magnetic field was increased, we 
observed periodic or quasi-periodic heterodyne voltage 
oscillation (“HV-oscillation”). The voltage changes re- 
producibly and periodically with magnetic field, up to 
the first critical field of the electrodes, Bd «1.6 Gauss. 
This is the critical field at which Abrikosov vortices begin 
to enter the electrodes. The period of the HV-oscillation, 
AB, can be defined as the distance between the adjacent 
principle maxima, as is illustrated in Fig. 2(a) by the hor¬ 
izontal arrow. Equivalently, the period can be defined as 
the separation between the so-called sweet spots. Some 
of these sweet spots, namely B 4 and B^ are indicated 
by the vertical arrows in Fig. 2(a). The sweet spots are 
the points equivalent to zero magnetic field. It should 
be reminded here that if the device is tuned to a sweet 
spot then it is insensitive, in the first order, to the flux 
noise, because dfoi/dB=0 [see Fig. 3(a)] and dVii/dB=0 
[see Fig. 2(a)]. 

At sufficiently low fields, when there are now vortices 
in the electrodes, the sweet spots occur periodically be¬ 
cause the critical current of the SQUID loop changes pe¬ 
riodically with magnetic field. In our case the design is 
such that the phase gradient created by vortices present 
in the electrodes couples to the SQUID loop. Therefore 
the exact periodicity of the sweet spots becomes broken 
when vortices begin to enter the electrodes at B > B^. 
Thus our device acts as a vortex detector. Entrapment 
of vortices inside the electrodes makes the transmission 
hysteretic with magnetic field up-down sweeps. An ex¬ 
ample of such hysteresis is clearly seen at the position 
marked “S” in Fig. 2(a). In addition to the period in¬ 
creases and the hysteresis, we also observe abrupt jumps 
in the transmission. One such jump is marked “J” in 
Fig. 2(a). The jumps indicate that the vortex enters in 
the near proximity of the qubit loop and thus should be 
strongly coupled to the qubit state. 

In our qubits, the effective Josephson energy is mod¬ 
ulated by the external magnetic field, but the period is 


set differently as compared to the typical SQUID-type 
device^^. Unlike in usual SQUIDs, the modulation of the 
critical current in the present case is driven mostly by 
the Meissner currents in the electrodes, and to a much 
lesser extent by the magnetic flux through the SQUID 
loop. This can be seen from the fact (See Table I) that 
the experimental period, AB, is much smaller than the 
period computed using the area of the superconducting 
loop, AByz- This is why the sensitivity of the Meissner 
qubit energy to the external field is higher compared to 
the usual split-j unction transmon^^. 

One can understand this new period by considering the 
phase constraints®’^®: 

61-62 + 2S{B) = 27rnv, (1) 

where 6 * 1,2 is the phase difference across each JJ, S{B) the 
phase difference generated in the thin-film electrodes by 
the Meissner current and defined as the phase difference 
between the entrance points of of the JJ bridges, i.e., 
between the bridges in which JJs are created. Here Uv 
is the vorticity. Since the current-phase relationship of 
JJs is single-valued, and the inductance of the JJ bridges 
and the electrodes is very small, the vorticity is always 
zero, riv = 0, in our SQUID-type devices, just like in 
common SQUIDs. The field-dependent phase accumu¬ 
lation is 6 {B)=f Vip{B) ■ dl, where W(p{B) is the phase 
gradient of the order parameter in the electrodes. The 
gradient originates from the Meissner (screening) current 
in the electrodes, if there are no vortices. At sufficiently 
high fields, at which vortices enter the electrodes, an ad¬ 
ditional contribution to the total phase gradient occurs 
due to the vortices. Following the Ref. 26, the magnetic 
period can be estimated as 



AB = 

[(*»)■' + 

f-Vl 

-1 

(2) 



\cXZj 

\YZj \ 


where 

the 

numerical 

coefficient 

c 

— 

(S/tt^) 

o(-l)”/(2n + 1)2 

« 0.74 can be 

found 


by solving appropriate boundary problem for the 
Laplace equation. We cautiously notice that the Eq. (2) 
may not be strictly applicable to our case because it 
was derived for mesoscopic electrodes, where A_l is 
much larger than X, the width of the electrode. In our 
samples A_l « 550 nm < A = 10 or 15 p,m. Meissner 
currents are stronger in the case of a relatively small 
perpendicular magnetic length, therefore this model may 
still provide a semi-quantitative estimate. To achieve 
a satisfactory agreement we will have to introduce 
corrections related to the field focusing effect. 

We emphasize that unlike a regular SQUID the period 
is not set by the area YZ enclosed by SQUID loop only, 
but rather by a much larger effective area YZ + cXZ; 
the dimensions re indicated in Fig. l(d)^®’^®. To com¬ 
pare those two effective areas, we plot in Fig. 2(b) the 
theoretical effective areas versus the “experimental” 
effective area A“ = $o/Ai? for five samples. Here AB is 
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the low-field, vortex-free period of the HV-oscillation for 
each sample. The black squares and the red triangles rep¬ 
resent, respectively, the theoretical effective area calcu¬ 
lated as =F Z (geometric SQUID loop area approach) 
and as Alfi=YZ + cXZ (Meissner current phase gradi¬ 
ents approach). The dashed curve represents the ideal 
case, The red triangles appear much closer to 

the ideal dashed line. Therefore the qubit energy is con¬ 
trolled mostly by the Meissner currents, which produce 
a strong phase bias of the SQUID loop. 

Another possible contribution to the observed ampli¬ 
fication of the magnetic sensitivity is the focusing of the 
magnetic field into the SQUID loop area by the supercon¬ 
ducting electrodes. Such focusing is also due to Meiss¬ 
ner effect, which is strong due the relatively small A_l. 
This field focusing effect enhances the magnetic field by 
a factor oi k = Bi/Bq > I. The ratio of the field i?i, 
enhanced by the field-focusing, to the applied field Bq is 
estimated in Appendix B. Thus the effective area set by 
the SQUID loop increases by a factor of k. To incorpo¬ 
rate the field focusing effect, we replaced YZ with kYZ 
in Eq. (2). The result is plotted in Fig. 2(b), showing an 
improved and now quite satisfactory agreement with the 
experimental results. 

Now, we turn to the magnetic field dependence of the 
qubit energy. For the spectroscopy, the qubit was excited 
with 2 /xs long saturation pulse, which was immediately 
followed by a few microsecond readout pulse. The exci¬ 
tation frequency was swept up—low to high frequency— 
with a fixed step size at a fixed magnetic field, and this 
process was repeated for equally spaced magnetic field. 
Fig. 3(a) shows the 2D color plot of the transmission as 
a function of the excitation frequency and external mag¬ 
netic fields. The color represents the heterodyne volt¬ 
age of the transmission of the cavity. The dashed line is 
a fit to the qubit spectrum with the following fit func¬ 
tion: /oi=/o\/| cos(7r(B - Bo)Arff/4>o)|, where /o, Bq 
and A^tt are the fitting parameters. We used the approx¬ 
imate relation for /oi=-\/8FjFc//i, where E,j=1llc/2e, 
and /c(B)=2/ci| cos(7r$/$o)|=2/ci| cos(7rBA,ff/$o)|. Bq 
and A^f! are offset field to account for the residual mag¬ 
netic field and the effective area, respectively. The best 
fit values of Bq and /q were —4.5 mG and 6.583 GHz. 
The best fit effective area was A^n = 55.2 /xm^, which is 
consistent with the = 55 /xm^ determined from the 
periodic oscillations of Fig. 2(a). 


B. Time domain measurement: low magnetic field 

Now we will look into the conventional time-domain 
measurements under the applied magnetic field. The 
time-domain measurements shown in Fig. 4 were per¬ 
formed to measure three time scales: relaxation time 
(Ti), phase coherence time (T 2 *) by Ramsey fringe, and 
phase coherence time (T 2 ) by Hahn spin echo. 

For the relaxation time measurement, we applied tt 
pulse (100 ns) first and then read out the qubit state af- 
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FIG. 3. (Color online) (a) Spectroscopy of Meissner transmon 
(Nl) as a function of applied magnetic field. This is raw data, 
(b) The parabola-like dashed line shows a phenomenological 
fit to the qubit transition frequency foi versus magnetic field 
B. 


ter the time interval At. For the Ramsey measurement, 
we applied two 7r/2 pulses (50 ns) separate by pulse sep¬ 
aration time Atft, and then readout was performed im¬ 
mediately after the second 7r/2 pulse. Similarly, in spin 
echo measurement, the measurement protocol was as the 
Ramsey protocol, except that an additional tt pulse was 
inserted right in the middle of the two 7r/2 pulses. The 
separation between the two 7r/2 pulses is denoted by Atg 
[see the inset in Fig. 4(c)]. 

Three time scales, Ti, T 2 and T 2 were ex¬ 
tracted by fitting data with exponential decay function 
exp(—t/to) for Ti and r 2 , and sine-damped function 
exp(—t/to) sin {2Tr+ (fo) for T|. In Fig. 5, we present 
/oi, Ti, T 2 and T 2 versus magnetic field measured in a 
vicinity of the sweet spot at B=0. The used weak mag¬ 
netic field, up to ~50 mG, is much weaker than the field 
needed to drive vortices into the electrodes. Thus here 
we discuss the vortex-free regime. 

We first examined the energy relaxation time, Ti, for 
both samples. They were measured separately in the 
same cavity which has its loaded lowest order mode 
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FIG. 4. (Color online) Time domain measurements of the 
N7 sample at i? = 7.5 mG. (a) Relaxation time measurement 
(Ti = 51 fj,s) (b) Ramsey fringe experiment (T)" = 18 fj,s) 
(c) Hahn spin echo experiment {T 2 — 27 ps). The red solid 
lines are the fits to the data. See the main text for the fitting 
functions. 


at ^8.42 GHz with a loaded (measured) quality factor 
Ql = 5000. As Fig. 5(a) and 5(b) show, at zero field the 
relaxation time Ti was substantially larger for the sam¬ 
ple N7 than for Nl. Furthermore, when a small magnetic 
field was applied, the energy relaxation time for Nl in¬ 
creased, while hardly any change was observed for N7. 
Both of these effects can be understood as consequences 
of the Purcell effect in which the rate of spontaneous 
emission is increased when the cavity mode to which the 
qubit couples lies close by in frequency. The excitation 
frequency of N7 (4.97 GHz) was further from the cav¬ 
ity frequency than was the excitation frequency of Nl 
(6.583 GHz) and this almost completely determines the 
difference in Ti. To see this we first compare the mea¬ 
sured ratio Fni/Fn? and the calculated ratio of Purcell 
rates for the two devices. (Since the next higher cavity 
resonance is more than 11 GHz above the fundamental, 
we ignore its contribution to the Purcell relaxation rate.) 
For the Purcell relaxation rate we use Fp = Ki{g/A)'^, 
where ki = lOcIQ'l is the cavity power decay rate, g is 
the qubit-cavity coupling rate, and A is the qubit-cavity 
detuning A = jwoi — wd = 27r|/oi — /c|. The ratio of 
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FIG. 5. (Golor online) (a) Magnetic field dependence of the 
qubit frequency (/oi), three measured time scales {T\,T 2 , 
and T 2 ) and two calculated time scales (Tp and Tf^*) at 
low magnetic field much smaller than the SQUID oscilla¬ 
tion period for the sample Nl. Tp (Purcell time) was cal¬ 
culated by Tp = l/Fp—inverse of Purcell rate, and Tf^* by 
1/T“i = 1/Tnp -F 1/Tp (see texts) (b) The qubit frequency 
and three measured time scales for N7. 


the Purcell rates depends only on the qubit-cavity fre¬ 
quency differences which are easy to measure. We find 
that the ratio of the qubits measured lifetimes for sam¬ 
ples Nl and N7 is Ti^ni/Ti^nt = 13 /rs/44 gs = 0.3, while 
the ratio of the calculated Purcell times, Tp = 1/Fp is 
(Ani/Ant)^ = 0.29. Since the measured and the esti¬ 
mated ratios are very close to each other, one can con¬ 
clude that the relaxation is Purcell limited. 

To confirm this conclusion we analyze the field depen¬ 
dence of the relaxation time. This is done using the 
qubit-cavity coupling, g — 130 MHz (see Appendix A) 
and the formula for the Purcell rate Fp = Ki(g/A)^. The 
increase in relaxation time for Nl with the applied B-field 
can be understood as a consequences of the Purcell effect. 
Indeed Pp ^ 1/A^ decreases with B because the qubit 
frequency decreases with increasing B-field thus causing 
the detuning A to increase. This makes a measurable dif¬ 
ference in Pp for Nl. For Nl we plot the relaxation time 
versus the B-field calculated as 1/Ti(i?) = Tp{B) -F Fnp 
[F ig. 5(a)]. Here Fnp = 27 r x 3.5 kHz is a constant Non- 
Purcell rate estimated at zero magnetic field. A good 
agreement between the data and the fit is observed, con¬ 
firming that the qubit is Purcell limited. Note that sam¬ 
ple N7 does not show a noticeable dependence on the 
B-field because its detuning value is high and therefore 
the Purcell effect is relatively weak. 

Finally we can estimate the internal relaxation rate of 
our devices assuming that the Purcell effect is eliminated 
by making the difference between the qubit frequency 
and the cavity difference sufficiently large. Such internal 
relaxation is represented by Fnp, as explained above. For 
Nl we obtain, at the sweet spot, Fp = 27r x 8.5 kHz and 
Fnp = 27r x 3.5 kHz, and, correspondingly Tnp = 45 /rs. 
For N7 we estimate, again at the zero field sweet spot, 
Fp = 27rx2.4 kHz and Pnp = 27rx 1.2 kHz, and, therefore 
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Tnp = 132 fj,s. This analysis reveals that the relaxation 
time could be above 100 /iS if it were not Purcell-limited, 
indicating that the coupling to energy absorbing defects 
in the circuit and qubit is low. 

We now consider the spin echo coherence time T 2 at 
B — 0 which measures the phase coherence attainable in 
a qubit process. T 2 is related to Ti by the constitutive re¬ 
lation l/r 2 = l/(2Ti) -I- 1/T,^, where is the dephasing 
time due to random fluctuations of the phase evolution 
rate of the qubit wave function. N1 had a much shorter 
T 2 than N7. Specifically, T 2 ^ni = 4.2 /rs and T 2 ,n 7 = 39 
/rs, almost ten times as large. Using the measured val¬ 
ues of Ti we obtain r,^^Ni = 3 p,s and = 70 /xs. 

We attribute the much longer dephasing time for N7 to 
the fact that the testing conditions were different. For 
measurements of N7, base-temperature copper powder 
filters were added to the input and output ports of the 
cavity. These are known to reduce stray-photon noise by 
providing attenuation at low temperature. Such photon 
noise may be responsible for the significantly lower de¬ 
phasing time seen in the N1 measurements. The photon 
noise can induce dephasing because of a strong ac-Stark 
shift^^. Meanwhile, the T 2 of sample N7 became shorter 
as magnetic field was applied. Such behavior can be ex¬ 
pected because an application of even a small U-field 
leads to a shift from the magnetic sweet spot. Thus the 
qubit becomes more susceptible to dephasing caused by 
flux noise. This effect is less visible in the sample N1 
because T 2 is already strongly suppressed by the stray- 
photon noise in this sample. 


C. Time domain measurement: high magnetic field 


Now we consider a different regime where sufficiently 
high magnetic fields creates vortices on the electrodes. 
We investigate the effect of Abrikosov vortices on the 
coherence times. Since the samples were zero-field-cooled 
there were no vortices in the electrodes to begin with. We 
managed to gradually increase the number of vortices by 
sweeping up the external perpendicular magnetic field. 
All measurements presented below have been made at the 
sweet spots, which occur periodically or approximately 
periodically with the magnetic field. 

In contrary to a single transmon, the Meissner trans- 
mons have advantage to allow us to detect the entrance 
of a vortex (or multiple vortices) into the electrodes. In 
actual measurements, we ramped up the magnetic field 
until the next sweet spot was reached [see Fig. 2(a)] and 
then performed the next series of time-domain measure¬ 
ments. Upon the event of vortex entrance, we observed 
two signatures: hysteresis of the transmission plotted ver¬ 
sus the magnetic field and a shift of the next sweet spot to 
higher magnetic field than would be expected if the pat¬ 
tern was exactly periodic. This happens because each 
vortex in the electrodes adds a phase gradient opposite 
to the one generated by the Meissner current [Fig. 2(a)]. 

In Fig. 6(a), we show the magnetic field dependence 


of three times scales and the qubit frequency for sample 
Nl. It is observed that Ti was enhanced from 10 p,s to 
14 /xs as the magnetic field was increased. The trend was 
observed up to about 2 Gauss, while at higher fields the 
trend was reversed. In the case of sample N7 [Fig. 5(b)], 
all three time scales stayed almost constant up to about 
2 Gauss and started to drop as the field was increasing 
further. This is explained by the fact that vortices be¬ 
gin to penetrate into the electrodes at the first critical 
field Bel ~2 G. They provide a radiation-free dissipation 
source and thus suppress the relaxation time significantly 
[see Fig. 6]. Meanwhile, the coherence time T 2 (and T 2 ) 
also decreased, mainly due to the reduction of Ti. Since 
the measurements were carried out at the sweet spots, 
the dephasing caused by non-zero dispersion of foi{B) 
was negligible. 

The entrance of vortices is confirmed by a compar¬ 
ison to a theoretical model. In Ref. 28 the first crit¬ 
ical field is approximated as i?ci=4>o/[27r^-v/2AXX], 
where Ax=2A(0)^/d represents the penetration depth of 
a thin film in perpendicular field. Using the relations, 
A(0 )=AlV1 -f and we estimate Bci=5.9 

G which is similar to the measured value, 2 G. Here 
^ is the coherence length, ^0 is the clean limit coher¬ 
ence length, A(0) is the bulk penetration depth, Al is 
the clean-limit penetration depth for Al. The parameter 
values are: ^0 = 1600 nm, I = 16.7 nm, Al = 16 nm, 
^=163 nm, A(0) = 158 nm. Ax = 552 nm, X = 10 /xm. 
The electronic mean free path I is calculated from the 
measured resistivity of the Al films forming the elec¬ 
trodes, pn = 2.4 X 10“®, according to Ref.^®, using 
lpn = 4:X 10-1® Um2. 

Our goal now is to achieve a quantitative characteri¬ 
zation of the non-radiative relaxation process caused by 
vortices. For this, we plot the relaxation rate F = 1/Ti 
versus magnetic field, B, in Fig. 7. The relaxation rate re¬ 
mains approximately constant at B < B^i and increases 
approximately linearly at H > B^i- The observed in- 


(a) 



Magnetic Field(G) 


(b) 



Magnetic Field (G) 


FIG. 6. (Color online) The qubit transition frequencies (/oi) 
and three times scales {Ti,T 2 , and T 2 ) were measured at the 
sweet spots over the wide range of magnetic field for the Nl 
(a) and N7 (b). 
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(b) 



Magnetic Field (G) 


FIG. 7. (Color online) The relaxation rates r=l/Ti versus 
magnetic field were plotted for both N1 (a) and N7 (b). All 
measurements have been done at the sweet spots. The linear 
fits (red solid lines) were shown over the magnetic field range 
where the vortices were present in the electrodes. 


crease of the relaxation rate per Gauss {dV/dB) was 78.5 
kHz/G for sample N1 and 43.7 kHz/G for sample N7. In 
what follows we suggest a model of non-radiative decay, 
which can explain these values. 

We suggest that the energy relaxation of the qubits 
is mainly due to the energy dissipation originating from 
the vortex viscous motion. The motion of vortices is 
initiated be the Lorentz force, which, in turn, is due to 
the currents generated by the qubit itself. The motion 
of vortices is overdamped due to the viscous drag force, 
which, per unit vortex length, is = —Cv, where ^ is a 
viscous drag coefficient for a vortex of unit length, and v 
is the vortex velocity^*^. This process is a non-radiative 
relaxation in which the qubit energy is dissipated as heat. 
We estimate this relaxation rate semi-classically, using 
the Bardeen-Stephen model^°. According to their model, 
the viscosity ( per unit length is C = 4’oBc2//On, where 
Bc 2 = 4>o/(27r^^), and ^ Here 3^2 is the second 

critical held of the thin-hlm A1 electrodes and is their 
normal-state resistivity. 

Let us estimate the energy relaxation rate, Tv, caused 
by one vortex via viscous damping. The rate of the en¬ 
ergy dissipation—dissipated power—is P = — (fv • v)d = 
where d is the thickness of him. Thus the energy 
relaxation rate of a transmon by a vortex can be evalu¬ 
ated as rv=P/(?iwoi) = (/v/C) d/hujQi, where Tiujqi is the 
energy stored in the hrst excited state of the qubit. The 
vortex is driven by the Lorentz force fL=J x <i)o (this 
is the force per unit length), where J is a supercurrent 
density. The supercurrent density is proportional to the 
total current, J = I/Xd, i.e., the current density magni¬ 
tude is approximated by the total current I divided by 
the cross section area of the electrode. A naive hrst guess 
might be that vortices in the electrodes should not move 
at all since the expectation value of the current generated 
by the qubit is zero, (J) = 0, both in the ground and in 
the exited state of the qubit. Yet we will see soon that 
the dissipation is proportional to (J^), which is greater 
than zero. 


The next step is to set /v = /l based on a reasonable 
assumption that the effective mass of the vortex and the 
pinning force are negligible. Gonsequently, we obtained 
the energy relaxation rate per vortex. 


(hujoi ’ 


( 3 ) 


Of course for the quantum states of the qubit the cur¬ 
rent and the current density are quantum variables which 
do not have definite values but should be viewed as quan¬ 
tum mechanical operators. The probability amplitudes of 
these quantities are defined by the wave function of the 
qubit. Thus under we understand the mean square 
of the current density, /{XdY, where the 

averaging is done for the excited quantum state of the 
qubit. Here / is the operator of the current in the qubit 
and X and d represent the width and the thickness of 
the electrodes. 

The model outlined above leads to the following esti¬ 
mates for the relaxation rate induced on the qubit by 
a single vortex: Tv = 89 kHz/vortex and Tv = 48 
kHz/vortex, for samples N1 and N7 correspondingly. The 
following set of sa mple- specific parameters has been used 
for sample Nl, sj(P) = 29.5 nA, J = 32.8 kA/m^, 
woi/27r = 6.583 GHz, C = 1.1 x IQ-^ N-s/m^, 3^2 = 12.2 
mT, pn = 2.4 X 10“® H-m, ^ = 163 nm, = 1600 nm, 
X = 10 pm, d = 90 nm, and I = 16.7 nm. For sample 
N7, all the parameters are the same, with the exception of 
= 18.9 nA, J = 21.0 kA/m^, and woi/27r = 4.972 

GHz. 

To compare the relaxation rate Tv, computed per a 
single vortex (see above), with the experimental relax¬ 
ation rates, dT/dB, measured “per Gauss”, we need to 
estimate the average rate of the vortex entrance, dN/dB. 
Then one can use a formula dT/dB — T^{dN/dB), which 
assumes that the relaxation rates of all vortices simply 
add up. 

A detailed analysis of various possibilities to estimate 
the vortex entrance rate are given in Appendix C. Here 
we briefly outline the most intuitive estimate. First, we 
define B^ as the sequence of magnetic helds correspond¬ 
ing to the sequence of the sweet spots, indexed by the 
integer n=0, 1,..., 27. Here 27 is the total number of the 
observed sweet spots [see Fig. 8(a)]. 

At low fields Bn (black open circle) increase linearly 
with n, due to exact periodicity of the HV-oscillation 
in the vortex-free regime. The linear fit [blue line in 
Fig. 8(a)] gives the value of the period, AB —0.2 G. For 
n > 7 the period becomes larger because vortices begin 
to penetrate. The new slope, and, correspondingly, the 
new period is AB + AB^ = 0.278 G, for sample N7. This 
best-fit value is obtained from the linear fit represented 
by the red solid line in Fig. 8(a). The enlargement of 
the period happens because vortices compensate, to some 
extent, the strength of the Meissner current. 

Now we calculate the difference between the consec¬ 
utive sweet spot fields — Bn-i- The result 













FIG. 8. (Color online)(a) The magnetic fields Bn at sweet 
spots (black open circle) are depicted as a function of n — 
index of the sweet spots, for sample N7. The blue and red 
solid lines represent a linear fit to the data for n < 6 and 
n > 7, respectively. The blue fitted line is extended for n > 7 
to show the expected Bn when no vortex penetration is as¬ 
sumed. (b) The difference AB„=B„ — Bn-i is shown on 
the left axis, while the number of vorticies ANn (defined in 
the text) provides the scale for the right axis, (c) The to¬ 
tal number of vorticies is calculated by summation of ANn- 
The dashed lines are periodically spaced. They provide a 
guide-to-the-eye to emphasize the stepwise characteristics of 
the increasing number of vortices in the electrodes. 

is plotted in Fig. 8(b). The result was then converted 
into some effective change of the vortex number, AN ^ 
{Bn — Bn-i) — AB, where AB is the distance between 
the sweet spots in the vortex-free low-held regime. Thus 
dehned AN should be considered as a function propor¬ 
tional to the number of vortices coupled to the qubit. 
But since the conversion factor is not well known, AN 
should not be considered as equal to the number of rele¬ 
vant vortices. Finally, we integrate AiV„ with respect to 
n. The result is shown in Fig. 8(c). This integrated func¬ 
tion exhibits clear steps, which we interpret as vortices 
entering the sensitivity area of the qubit. The steps are 
made more noticeable by placing the horizontal dashed 
lines. The spacing between the lines is made constant 
and they serve as guide to the eye. The total number of 
effectively coupled vortices equals the number of steps, 
i.e. equals 6. These 6 vortices have entered over the in¬ 
terval of 5.6 Gauss. Thus the effective entrance rate is 
estimated as dN/dB « 1.07 vortex/G. 

Remember that the experimental relaxation rates per 
Gauss dr/dB, obtained from Fig. 7, are 78.5 and 43.7 
kHz/G for samples N1 and N7. Now these values need 
to be divided by 1.07 vortex/G, which is the rate of the 
vortex entrance. Thus we conclude that the experimen¬ 
tal relaxation rates are Fv = 73 kHz/vortex and Fv = 41 
kHz/vortex, for samples N1 and N7, calculated using 


Fv = {dr/dB)/{dN/dB). These values are in excellent 
agreement with the theoretical estimates 89 kHz and 48 
kHz/vortex. 

To understand the signihcance of the obtained results 
for the macroscopic quantum physics, it is instructive 
to compare the relaxation rates generated by individual 
vortices to the theory of localization caused by dissipa¬ 
tive environment. According to the Caldeira and Leggett 
(CL) theory^^, the particle wave function should be ex¬ 
ponentially localized with the localization length scale 
estimated as — h/r], where r] is the viscous drag co¬ 
efficient of the macroscopic quantum particle coupled to 
the environment. This theory was later generalized to 
the case of periodic potentials in Refs.^^’^^. The conclu¬ 
sion of these papers was that if the period of the potential 
is larger than xcl then the particle becomes localized in 
one of the wells in the limit of zero temperature. On the 
other hand, if the period is smaller than xql then the 
particle can tunnel from one minimum to the next one 
even at zero temperature. It is important that the scale 
of the environmental localization, xcl, is independent of 
the amplitude of the modulation of the potential energy. 
Thus it can provide a useful estimate even if the potential 
is approximately flat, which is the case for the Abrikosov 
vortex in the A1 film in our devices. The viscous drag 
coefficient for a single vortex is rj = dC, = 9.9 x 10“^^ 
N-s/m, where d = 90 nm is the vortex length, which is 
set by the film thickness d. Then the Caldeira-Leggett 
(CL) localization scale is xcl = y/h/rj = 2.6 x 10“® m = 
2.6 nm. This is the scale up to which the wave function 
of a vortex can spread. If the spread is larger than xcl 
then the coupling to the environment causes the wave 
function collapse. 

The CL localization scale should be compared to the 
estimated smearing of the wave function of the center- 
of-mass of the vortex, generated by fluctuations of the 
supercurrent and the corresponding fluctuations of the 
Lorentz force. The smearing of the wave function can be 
estimated as follows. The root-mean-square (rms) value 
of the Lorentz force is Fl = d • /l = 6.1 x 10“^® N. (Here 
we consider the example of sample Nl. All estimates for 
sample N7 are very similar.). Then, assuming viscous 
motion, the rms velocity is u = Fi^/rj = 6.2 x 10“^ m/s. 
Therefore, the rms quantum fluctuations of the vortex 
center position, Xv, can be estimated as Xv = v/uiqi. This 
relation would be exact if the motion of the vortex would 
be described by a classical trajectory, in response to a 
harmonic drive, such that the deviations from its point 
of equilibrium would be proportional to sin(a;oit). In the 
case considered the vortex is not described by a classical 
trajectory, since it should act as a quantum particle at 
time scales shorter than the relaxation time of the qubit 
and also because the force is generated by a quantum 
superposition of currents with opposite polarities. But 
we assume that the relationship between the quantum 
fluctuation of the position and the velocity is approxi¬ 
mately the same as in the case of a classical harmonic 
motion. Such assumption is motivated by the natural 
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expectation that the vortex should behave as a damped 
quasi-classical particle. Thus we can now evaluate the 
rms smearing, of the wave function of the vortex cen¬ 
ter induced by the quantum fluctuations of the Lorentz 
force. The result is Xy = (Tl/? 7 )/woi = 1.5 x 10“^^ 
m = 1.5 pm. Here we have used the qubit frequency 
/oi = a;oi/27r = 6.58 GHz for sample Nl. The conclusion 
is that the quantum uncertainty of the vortex position is 
much smaller than the CL localization scale, <C xcl, 
namely xc'l/x^, = 1700. Such uncertainty is achieved 
within one period, Tgi = 27r/woi- 

The CL scale provide the maximum value of the rms 
smearing of the wave function. When such level of smear¬ 
ing is achieved, the wave function collapses and the 
qubit experiences a decoherence event. The number of 
complete phase revolutions of the qubit, iVdch, which is 
needed to achieve the CL scale, at which the probabil¬ 
ity of decoherence becomes of the order of unity, can be 
estimated assuming that the wave function of the vor¬ 
tex spreads similar to a diffusive random walk. Then 
IVdch = = 3.0 x 10®. Thus the corresponding 

decoherence time can be estimated as tdch -VdchVoi- 
This heuristic argument can be made more precise us¬ 
ing Eq. (3). From that equation, using the relations 
listed above, one obtains tdch = (l/27r)^A^dchToi = 12 /xs. 
Finally, the estimate for the relaxation rate, added to 
the system due to one vortex coupled to the qubit, is 
Tv = 1/tdch = 83 kHz. This analysis leads to a con¬ 
clusion that the relaxation rate induced by one vortex 
can be linked to the ratio of the smearing of its wave 
function and the Caldeira-Leggett dissipative localization 
scale. The above estimate shows that the wave function 
collapse occurs typically at the time when the wave func¬ 
tion smearing reaches the Caldeira-Leggett localization 
scale. 


III. SUMMARY 

We demonstrate the operation of Meissner qubit, 
which is a transmon qubit strongly coupled to the Meiss¬ 
ner current in the electrodes. The periodicity in magnetic 
field was set mainly by the width of the electrodes rather 
than the SQUID loop area. Both the frequency and 
the time domain measurements were performed using the 
circuit-QED architecture. The three time scales of Ti, T| 
and T 2 were measured as functions of the applied mag¬ 
netic field. The increase of the relaxation rates was at¬ 
tributed to the radiation-free dissipation associated with 
the viscous motion of Abrikosov vortices pushed by fluc¬ 
tuations of the Lorentz force. Each vortex can exist in 
a quantum superposition of different position states for 
^10-100 /xs, but eventually causes the wave function of 
the coupled qubit to collapse. The collapse happens when 
the smearing of the vortex center becomes of the order of 
the Caldeira-Leggett dissipative localization scale. The 
presented Meissner qubit provides an effective and con¬ 
trolled coupling of the qubit to Abrikosov vortices. Such 


coupling provides a new tool to study vortices, which 
can eventually be applied to vortices harboring Majorana 
states. 


ACKNOWLEDGMENTS 

This work was supported by the National Science 
Foundation under the Grant No. ECCS-1408558 (A.B.) 
and ECCS-1407875 (A.L.). 

Appendix A: Device and Measurement 

The Meissner qubits were fabricated on a c-plane sap¬ 
phire using a novel modification of the double-angle evap¬ 
oration technique, now achieved in ultra-high vacuum in 
a molecular beam epitaxy (MBE) growth system. The 
device design is shown in Fig. 1(d). The patterns were 
defined by electron beam lithography on the bi-layer of 
MMA EL-13 and ZEP 520 A7 in eLine Raith system, 
and after development the exposed surface was cleaned 
to remove MMA residue by both dry and wet etching— 
oxygen plasma by RIF (reactive-ion etching) and BOE 
(buffered oxide etch). The first and second layers of alu¬ 
minum films, each with 45 nm thickness were deposited 
with background base pressure of 10“^^ Torr. The ox¬ 
ide layer was formed by an exposure to Ar /02 mixture 
(10% O 2 ) under the proper conditions of pressure and 
time calibrated for critical current density of the JJs. 
Each large rectangle [marked A1 and A2 in Fig. 1(b)] 
acts as a radio-frequency (RE) antenna and has dimen¬ 
sion of 250 X 500 /xm^. The spacing between the nearest 
edges of the antennas is 25 /xm [Fig. 1(b)]. The antennas 
are bridged by two A1 thin-film rectangles [called “elec¬ 
trodes” and marked El and E2 in Fig. 1(c)] and two JJs, 
connecting the electrodes and forming a SQUID-like loop 
[Figs. 1(c),1(d)]. 

We will discuss two representative devices denoted by 
Nl and N7. The qubit transition frequency /oi of Nl 
at zero magnetic field was 6.583 GHz, the Josephson en¬ 
ergy of both junctions taken together was U™“=19.4 
GHz, the corresponding net critical current was 41 nA, 
and the Coulomb charging energy—mostly associated 
with the electric capacitance between the antennas— 
was Uc=0.307 GHz. Thus the ratio of the two en¬ 
ergy scales was EJ^^^/Ec=63.2 for device Nl. For the 
qubit N7 the same type of parameters were /oi =4.970 
GHz, E^^^=ll.l GHz, the corresponding critical cur¬ 
rent at zero field was 23.4 nA, Fc=0.318 GHz, and 
Ef^^ / Ec=M.9. 

Each device was mounted in a 3D rectangular copper 
cavity with the dimension 28 X 22 X 5 mm® [Fig. 1(a)]. The 
qubits were coupled to electromagnetic cavity modes and 
in turn the cavity was used for qubit readout. The sam¬ 
ples were positioned in the center of the cavity where the 
dipole coupling of the fundamental mode to the qubit 
is maximized, i.e., at the electric field anti-node. The 
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TEoii mode of the empty cavity occurred at 8.679 GHz 
with the internal (unloaded) quality factor Qi=8000, ex¬ 
ternal quality factor Qe=14000, and loaded (’’experimen¬ 
tal”) quality factor (5 l= 5000^^. An asymmetric coupling 
(larger at output) was used to maximize the signal to 
noise ratio^®. The bare cavity mode at TEqh for 
sample N1(N7) was 8.403 GHz (8.435 GHz), and the 
coupling strength was approximately 130 MHz^^ 

for both samples. 

The cavity was mounted in a ^He/^He dilution refrig¬ 
erator (S.H.E. Corp.) with the base temperature of 45 
mK. The cavity was enclosed in a cylindrical aluminum 
Faraday cage, whose inner walls were coated with black 
infrared-absorbing material^®. The A1 cage is intended 
to prevent external stray photons from reaching the sam¬ 
ple. The A1 cage helps also to protect the sample from 
the influence of external stray magnetic helds, thanks 
to the Meissner effect. A cylindrical cryogenic ^-metal 
(Amuneal) cylinder was placed concentrically around the 
A1 cage for additional magnetic shielding. The magnetic 
field was applied perpendicular to the substrate from an 
external home-made superconducting solenoid attached 
at the bottom of the copper cavity, i.e. inside the A1 cage 
and the /r-metal shield. 

For microwave transmission measurement, the input 
and output transmission lines were connected in series 
with a chain of cryogenic microwave components, includ¬ 
ing attenuators, isolators (PAMTEQ), a commercial low- 
noise HEMT (high-electron-mobility transistor) amplifier 
(Low Noise Factory, LNF-LNC6-20A) and low-pass Hl- 
ters. For noise filtering, commercial low-pass filters (K&L 
Microwave, 6L250-12000) and home-made stainless steel 
powder filters (3 dB at cutoff « 8 GHz) were inserted in 
both input and output of the copper cavity. 

The measurements were performed using the circuit 
QED technique in high-power regime^°’^^. For the spec¬ 
troscopy and qubit readout, RF square pulses were cre¬ 
ated, added together and fed into the cavity. The pulses 
were shaped by mixing a continuous microwave tone 
(Agilent E8267C or HP 8341B) and a square voltage 
pulse from an arbitrary waveform generator (Tektronix 
AWG520) using pairs of rf mixers (Markis, M8-0420). 
The transmitted readout signal was down-converted to 
25 MHz intermediate frequency (IF) signal by heterodyne 
demodulation, and the IF signal was acquired to read the 
amplitude by a high-speed digitizer (Agilent, U1082A- 
001). For qubit state readout, we adjusted the power of 
readout pulse (a few gs long) to maximize the contrast 
in transmitted microwave amplitude for the ground and 
first excited states. 


Appendix B: Calculation of ft: Magnetic field 
focusing effect 

In this appendix, we will show how to calculate ft. Con¬ 
sider a two-dimensional array of square superconducting 
films with magnetic held Bq applied perpendicular to the 



FIG. 9. (Color online) An array of square superconducting 
aluminum films. A magnetic field is applied perpendicular to 
the surface of the films. The enclosed area by dashed line 
indicates one unit cell. 

in-plane of the array as shown in Fig. 9. The squares 
represent the electrodes in actual Meissner qubits. We 
consider a unit cell enclosed by a red dashed line to cal¬ 
culate ft. When the magnetic held Bq is applied over the 
unit cell, the magnetic held inside the hlms is expelled 
by Meissner effect, and thus the magnetic held (i?i) in 
the hatched area is enhanced by a factor of ft = B^/Bq. 
ft is calculated in the following way. We denote Aq to 
the area of a unit cell and Ai the hatched area. The 
magnetic hux in one unit cell is 4) = BqAq = BiAi, so 

B^ Aq (X + Y)^ {X + Yf 

Bo Ai (x + vy-x^ V(2X + Y) ^ ’ 

We note that ft > 1, i.e., held focusing effect. 

Appendix C: Analysis of the average number of 
vortices entering the electrodes per one Gauss of the 
external field, dN/dB. 

It is not possible to determine dN/dB exactly. There¬ 
fore we outline three diherent approaches below. The 
hrst two methods, which are very similar, will provide 
the higher bound for dN/dB, while the third one will 
give the lower bound for dN/dB. 

First, we dehne B^ as the sequence of magnetic held 
values corresponding to consecutive sweet spots, indexed 
by the integer n=0, 1,..., 27. Here 27 is the maximum 
number of the sweet spots measured on sample N7. The 
sweet spots correspond to the minima of the HV-signal 
shown in Fig. 2(a). 

At low helds [for n < 6 in Fig. 8(a)], (black open 
circle) increase linearly with n, as is expected for the sit¬ 
uation in which the sweet spots occur periodically with 
magnetic held. Such exact periodicity is observed only 
in the low held regime, when there are no vortices in the 
electrodes, i.e., < Bd for n < 6. Since the oscillation 

is perfectly periodic in this vortex-free regime, the posi¬ 
tions of the sweet spots of the qubit can be approximated 
as Bn = AB-n+Bo, where AB is the unperturbed period 
of the HV-oscillation and Bq is the position of the zero’s 
sweet spot. The linear ht [blue line in Fig. 8(a)] in the 
low-hied regime gives the value of the period, AB =0.2 
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G for n < 6. Note that in such representation the best 
linear fit provides the averaged period, AB = dBn/dn. 

In what follows we discuss the regime occurring above 
the critical field, with vortices entering the electrodes as 
the magnetic field is swept up. Since the slope of the 
Bn versus n dependence changes significantly at n = 7, 
therefore B-j « B^i. For n> 7 the period is larger com¬ 
pared to the unperturbed case with no vortices. This 
is because the current of each vortex is opposite to the 
Meissner screening current. Thus the total phase bias im¬ 
posed on the SQUID loop by the electrodes^® increases 
slower with the magnetic field if the number of vortices 
in the electrodes increases with magnetic field. Thus the 
sweet spots tend to occur at higher field values. The 
new dependence of the position of the sweet spots versus 
their consecutive number is still approximately linear [see 
Fig. 8(a)], but the slope is larger. The formula for the 
sweet spot sequence becomes « {AB + AB^)n -\- B'^, 
where AB^ is the value by which the average interval be¬ 
tween the sweet spots is increased, due to the continuous 
increase of the number of vortices N in the electrodes. 
The new slope, and, correspondingly, the new period is 
AB + ABv = 0.278 G, in the present example of sam¬ 
ple N7. This best-fit value is obtained from the linear fit 
represented by the red solid line in Fig. 8(a). 

The enlargement of the period ABy — 0.078 G is at¬ 
tributed to the additional phase bias®® induced on the 
SQUID loop by the vortices entering in the electrodes as 
the field is swept from one sweet spot to the next one. 
This change in the period can now be used to estimate 
the number of vortices entering the electrodes within one 
period of the HV-oscillation. Let us Atpv be the average 
phase difference imposed by a vortex on the SQUID loop 
and AN be the average number of vortices entering the 
electrodes during each period. Then the additional phase 
difference Atpvv accumulated during each period due to 
all newly entered vortices is A(pvv = A(pv • AN. Since the 
vortex current opposes the Meissner current, the phase 
difference Aipw is opposite in sign to the phase A<pm 
imposed by the Meissner current. Therefore the total 
phase bias generated within one period can be written 
as Atpt = A(^m — At^v As in any SQUID-based device, 
one period of oscillation corresponds to the total phase 
change by 27r. Thus one has to require Aip^ = 27r and so 
A(/5m = 27r -|- Aip^AN. 

The value of A(pM is determined by the phase-bias 
function 26{B), which is the function defining how much 
phase bias is produced by the applied magnetic field 
B [see Eq. (1)], taking into account only the Meissner- 
current-generated phase gradients in the electrodes®®’®®. 
Since in the vortex-free regime the period equals AB and 
the phase bias should change by 27r to complete one pe¬ 
riod, and also because SB is a linear function of B, we 
can write 

2S{B) = 2tt{B/AB) (CI) 

for the vortex-free regime. As vortices begin to en¬ 
ter the electrodes, the period increases, on average, to 


AB -I- Ai?v, as was discussed above. Again, the total 
phase generated within one period should be 27r. Thus, 
when vortices are entering we can write 26{AB + ABy) = 
27r -I- AcpyAN. The last term is needed since at the end 
of one period the phase generated by the Meissner cur¬ 
rent has to be larger than 27r by as much as A(pvAN, 
to compensate the opposite phase bias, Aip^AN, gen¬ 
erated by the newly entered vortices. Remember that 
the function d{B) is defined by Eq. (Gl). Therefore 
2tt{AB -\- ABy)/AB = 27r -I- Aip^AN. Finally we get 
the formula AN = {2tt/ A‘.p^){AB^/AB)^ which defines 
the average number of vortices entering the electrodes 
per one period of the HV-oscillation. Thus, for exam¬ 
ple of sample N7, the number of entering vortices, per 
one period, is AN = 2. In this estimate we have used 
ABv = 0.078 G, AB = 0.2 G, and Aip^ = 1.22 (to be 
discussed in the following paragraph). Now we are ready 
to make an estimate of the average number of vortices 
entering the electrodes as the applied field is changed by 
one Gauss, which is dN/dB = AN/{AB + AB^) = 7.2, 
assuming that N{B), the number of vortices versus mag¬ 
netic field, is linear. 

The estimate above required us to make an assump¬ 
tion that each vortex generates a phase bias Aip^ = 1.22, 
on average. To justify this, consider one vortex in 
the center of one of one of the electrode, illustrated 
Fig. 1(d). Then, according to Ref.®®, A(^v is equal 
to the polar angle ©v subtended by the line connect¬ 
ing the entrance points of the two bridges leading to 
the JJs forming the SQUID. Therefore, for our samples, 
A(fy = ©v = 2 tan“^[2Z/(25 — V)] [Y and Z are shown 
in Fig. 1(d)]. Therefore, for Sample N7, the phase bias 
imposed by one vortex located in the center of the elec¬ 
trode is Aipv = 1.22 rad. The geometry of the electrodes 
of the sample N1 is very similar to N7, therefore we will 
assume the same phase bias per vortex for the sample 
Nl. 

Of course, the number of vortices entering the elec¬ 
trodes exhibits fluctuations. In what follows we present 
a different approach to analyze and average out these 
fluctuations. First, we calculate the difference between 
the consecutive sweet spot fields Ai?„=R„ — Bn-i- 
The result is plotted in Fig. 8(b). With these nota¬ 
tions ABn — AB is the increase of the measured pe¬ 
riod above the vortex-free period AB. The number 
of vortices entering the electrodes between two adja¬ 
cent sweet spots, AiV„, can be estimated using formula 
ANn/{ABn — AB)=AN/ABy because AN =< AA„ > 
and Ai?v =< Ai?„ — AB >. Here it is assumed that the 
number of vortices entering per period is linearly pro¬ 
portional to the period. Thus computed number AN is 
shown on the right axis in Fig. 8(b). Finally, we inte¬ 
grate ANn with respect to n to output the total number 
of vortices, A„, versus the sweet spot index n. The re¬ 
sult is shown in Fig. 8(c). From this plot one can esti¬ 
mate that that 38 vortices enter as the index is increased 
by 20 (remember that An = 1 corresponds to one pe¬ 
riod). Thus one obtains 1.9 vortices per period. Since 
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the period equals 0.278 G, on average, one estimates that 
dN/dB =6.8 vortex/G, in the case of sample N7. Since 
N1 has about the same size of the electrodes as N7, we 
use the same conversion factor dN/dB for Nl. 

Now we discuss our third approach to estimate the 
vortex entrance rate dN/dB — (iV„+i — Nn)/NB. This 
approach is based on the observation that the function 
computed by the algorithm outlined above, exhibits 
a stepwise increase as the magnetic field is swept linearly 
[Fig. 8(c)]. The steps are made more noticeable by plac¬ 
ing the horizontal dashed lines. The spacing between the 
lines is constant and they serve as guide to the eye. The 
step size turns out almost constant. We speculate that 
each step corresponds to the entrance of a single vortex 
which is effectively coupled to the qubit. This scenario 
assumes that not all vortices present in the loop are suf¬ 
ficiently well coupled to the supercurrent generated by 
the qubit but only those which enter the area near the 
SQUID loop. The physical reason for this is the fact that 
the current tends to be concentrated near the edges due 
to the Meissner effect. At the same time, many other 
vortices get pushed in the middle of the electrode, thus 
making their impact on the qubit very minimal. It is 
naturally expected that vortices entering the electrodes 
near the loop would make a relatively large impact on the 
change of the period, Bn+i — B^ and therefore can cause 
an sharp increase in the estimated Nn number. Thus the 
steps apparent in Fig. 8(c) represent the vortices effec¬ 


tively coupled to the qubit, and only these vortices are 
relevant for our estimate of the relaxation rate. In this 
scenario, the total number of effectively coupled vortices 
equals the number of steps, i.e. equals 6. These 6 vortices 
have entered over the interval of 5.6 Gauss. Thus the ef¬ 
fective entrance rate can be taken as dN/dB=6/5.6=1.07 
vortex/G. This estimate provides the lower bound for 
dN/dB. 

Remember that the experimental relaxation rates per 
Gauss dV/dB, obtained from Fig. 7, are 78.5 and 43.7 
kHz/G for samples Nl and N7. Now these values need 
to be divided by 6.8, which is the rate of the vortex 
entrance, dN/dB. Thus we conclude that the exper¬ 
imental relaxation rate is Fv = 11.5 and Fv = 6.4 
kHz/vortex, for samples Nl and N7, calculated using 
Fv = {dV/dB)/[dN/dB). These values are somewhat 
smaller than the theoretical estimates of 89 kHz and 48 
kHz/vortex. This fact serves as indirect evidence that 
the number of vortices is overestimated. 

With the conversion factor obtained by our third 
method, based on the observation of the steps, the ex¬ 
perimental relaxation rates per vortex become 73 and 41 
kHz/vortex for sample Nl and N7. These values are in 
good agreement with the calculated values, 89 and 48 
kHz/vortex, for samples Nl and N7. Thus the approach 
based on the step counting (i.e., our third method) ap¬ 
pears to be the most accurate for the estimation of the ef¬ 
fective number of vortices effectively coupled to the qubit. 
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